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l ■ Abstract. The N-dimensional generalization of Bertrand spaces as families of Maximally 

superintegrable systems on spaces with nonconstant curvature is analyzed. Considering the 
classification of two dimensional radial systems admitting 3 constants of the motion at most 
r^ ' quadratic in the momenta, we will be able to generate a new class of spherically symmetric 

Oh| M.S. systems by using a technique based on coalgebra. The 3-dimensional realization of 

these systems provides the entire classification of classical spherically symmetric M.S. systems 
admitting periodic trajectories. We show that in dimension N > 2 these systems (classical 
Cy ^ and quantum) admit, in general, higher order constants of motion and turn out to be exactly 

solvable. Furthermore it is possible to obtain non radial M.S. systems by introducing projection 
of the original radial system to a suitable lower dimensional space. 
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Q"^ I 1. Introduction 

'nI" ■ A Maximally Superintegrable (MS) system in classical mechanics is an integrable N-dimensional 

■^ \ Hamiltonian system which is endowed with the maximum possible number of 2N — 1 functionally 

^^ ' independent integrals of motion. In the last decades there has been a growing attention towards 

this class of systems due to their physically relevant properties both in classical and quantum 
mechanics : they are conjectured to be exactly solvable [3] , they are often multiseparable (the 
system can be separated in more than one coordinate system), the (M.S.) classical systems 
^ , exhibit periodic trajectories for all of their bounded motion and the quantum mechanical ones 

have degenerate spectrum for bound states. Because of these remarkable properties we find 
(M.S.) systems as successful models in many areas of physics such as in condensed matter 
physics as well as atomic, molecular and nuclear physics see e.g. [1] [5] [6] and reference therein. 
Superintegrable systems are very rare, in the sense that, under certain assumptions it is possible 
to classify completely all the possible systems in the (M.S.) family : The most celebrated example 
dates back to the 19th century and is a consequence of the Bertrand's theorem P]. Given a 
Hamiltonian system defined on a flat space, then the only central potentials admitting the (M.S.) 
property are 

H=pl + ?^ + V{r), y(r) = {^,^V}, 
the harmonic oscillator or Kepler-Coulomb potential. These two systems are also characterized 



by having quadratic constants of the motion /'^^ 



N 



{H,I} = : /(r,0,p„p,)W = J]^a(r,0)„,,pry„ N = 2. 

ra=0 i=0 

Over the last century many other classifications have been made, above all considering 
superintegrability characterized by constants of the motion I^ ' of a fixed order N defined 
on both Euclidean and non Euclidean spaces. The literature about it is simply too vast to 
cite all the important contributions, however a quite exhaustive list can be found in the review 
paper [J] where superintegrable systems with constants of the motion up to the third order are 
considered. The main goal of the paper is analyzing the classification of all the M.S. systems on 
general N-dimensional Non-Euclidean spaces characterized by radial symmetry and admitting 
bound states, which, in the 3-dimensional case coincide with the so called Perlick systems [2j. 
This class of systems has already been introduced in a series of papers by Ballesteros et al. and 
the author [8] [9] [10] |11] in which the Maximal Superintegrability of all the radial systems 
defined on Perlick's space-times is proven in their classical version and for few of them also 
in the quantum one. In this paper we will present a strategy to obtain systematically all the 
constants of the motion /'^' (which in general are of higher order A^) for these systems in 
spaces of arbitrary dimension both in classical and quantum mechanics. In doing this we will 
emphasize the special role played by the systems admitting quadratic constants of the motion 
H, j'-^': The classification of quadratic (M.S.) systems with radial symmetry on two dimensional 
non-Euclidean spaces admitting periodic trajectories can be expressed in terms of two families 
which can be regarded as deformations of the Kepler and Harmonic oscillator systems |12] |13j . 
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2(1 + A2r4 - 25r2) ^^"^ ' r^' ' 1 + A^r^ - 25r2 ' 
We will show how the canonical transformation 

r = r'P;Pr = -^Pr'] 

9 = f3e;pQ = —- 

/5 



introduces a new angular parameter f3 which, after a proper embedding procedure into a 3 



dimensional space Pq ^ Pg + . t^ > turns the 2-D quadratic systems into 



^I = 0R2 -^Pr' + -I2+ ,2 2 J + ^i"-^ - ^^^ ' 

2p"^ r'^ r'^ sm 6' 

/2(/-2/^ - A2r-2/3)2 ^ p2^ p2^, ^ 

^11 2(r'-2/3 + A2r'2/3 - 25) ^^^' r'2 + r'2 sin2 6' ^ r'-2/3 + A2r2/3 - 26 ' 

This is the entire Perlick family, which in contrast to the two dimensional case, is characterized 
by having in general higher order constants of the motion /'^' . As we will show in section @ 
this fact can be understood in terms of the sl2 coalgebra symmetry jl4| which is common to 
all the Perlick Hamitlonians regardless the dimension of the space in which they are embedded. 
This is in contrast with the representation of its algebra generators, which generally depends on 
the dimension of the space we are considering and that, in the particular case of a 2-dimensional 
space, admits always a quadratic representation for all the first integrals. 



Moreover the identification of tlie coalgebra structure for tliese systems opens a way to obtain a 
wider class of higher order superintegrable systems by considering just different representations of 
the same algebra. As a consequence of this we will show explicitly how it is possible, by choosing 
different representations of the coalgebra elements, to transform the family of Perlick systems 
into the family of the well known TTW systems and its extensions [15] |16] |17] . These exhibit 
higher order constants of the motion on the same footing of the Perlick systems. To conclude 
the paper we will discuss the quantum version of all these systems stressing the role of the SUSY 
QM (Super Symmetric Quantum Mechanics) that can be fruitfully used to get the solution and 
the superintegrable quantization of the MS Hamiltonians and its first integrals. The paper is 
organized as follows: In Section 2 we will give a survey of radial superintegrable systems on 
two dimensional non-Euclidean spaces. In Section 3 we will introduce the 3-dimensional Perlick 
systems stressing their connections with the classifications of the quadratic two dimensional 
superintegrable systems introduced in the Section 2. In Section 4 we will analyze the integrability 
properties of the Perlick systems in terms of the coalgebra sl2. In Section 5 we will extend the 
coalgebra approach to obtain non-radial higher order M.S. systems. Finally in Section 6 and 7 
we will provide the quantum version of all the systems previously analyzed in sections (1-5). 

2. Two dimensional Maximal Superintegrable radial systems 

It is very well known that for any classical N-dimensional radial Hamiltonian system all the 
trajectories describe a planar motion. Their quantum analogues have a spectrum which depends 
on just two quantum numbers. This fact can be explained in terms of superintegrablility claiming 
that every N-dimensional radial system has a " coalgebra-symmetry" |14j which makes the system 
itself quasi-maximally superintegrable, namely it has at least 2N — 2 integrals of motion whose 
main effect is that of reducing its dynamics to that of a 2 dimensional system. In the light 
of these considerations let us start our analysis giving a brief review of the classification of 
two dimensional M.S. systems. In particular we will show how the entire family of radial M.S. 
systems admitting periodic trajectories can be obtained from a special class of M.S. systems 
which is characterized by the fact of having at most quadratic integrals of the motion. Let us 
introduce a general two-dimensional Riemannian space whose metric is given by: 

ds'^ = gijdx'^dx^, i,j = 1,2 (1) 

The Hamiltonian describing a particle moving on such a space and subjected to a potential 
V{x^,x'^) has the following form for a classical system: 

H=-g''^PiPj+fiVixi,X2), (2) 

while for a quantum mechanical system we have: 

H = --L_o.(^^g^,JOj) + fiV{xuX2). (3) 

The form of the Hamiltonians introduced above induces a preliminary discussion on what are 
the possible Riemannian spaces admitting a M.S. Hamiltonian system. If we exclude from this 
discussion the trivial ones like the Euclidean space and the constant curvature space then the 
analysis reduces to the classification of all those Riemannian spaces with a non constant scalar 
curvature whose associated Hamiltonian admits three functionally independent constants of the 
motion. A significant result in this sense was obtained by Koenigs [IS] in 1889. He classified 
all those Riemannian spaces whose Hamiltonian admits three constants of the motion / at most 
quadratic in the momentum. 



More explicitly such a condition in classical and quantum mechanics is respectively equivalent 
to: 

I = a{xu XjY^^p^pj + /(xi, X2) : {H, /} = (4) 

i = a{xi,Xjy'^d,dj + b{xi,X2ydi + f{xi,X2): [hJ] = 0. (5) 

Where the curly bracket {,} stands for the Poisson brackets and the square one [, ] for the Lie 
commutator. The list of spaces which satisfy at these constraints are called Darboux spaces. 
This list contains 4 spaces which we report here in diagonal form as introduced in |20] |21] : 



(I) 2x{dx'^ + dy^) (6) 

(11) ^^{dx^ + dy^) (7) 

(///) ^^^{dx' + dy') (8) 

/T^.x 2acosh(2x) + 6 , , 9,9. 

These metric spaces can be recast in a radial conformally flat reference frame by defining the 
transformation {x,y) — )• (r = ^,y = 9), e = ±1. The Riemannian metric trasforms into 
fix){dx^ + dy2) ^ /(£W) (^^2 + ^2^^2^ 

(I) ^iM#i!:l(dr^+rW) (10) 

{III) "-^^^(dr^+rV) (12) 



^2^2e^2e _ j.-2e-\2 



m - '..ZT'-Zr id-' + r^de^)- (13) 



Let us note that the Darboux III and Darboux IV can be described as particular cases of a more 
general Riemannian space : 



9 2c(r-2^ + AV^-2<5),_^9 , 9,,2^ 



ds' = Uir-^e _ X2r2e^2 ' (dr^ + ^'d^ (14) 



It is straightforward to see that the Darboux IV metric coincides with ()14p if we set the 

parameters equal to c = — 2^^^ = "20"^^ = -^1- The previous change of parameters 
is not defined for A = 0. This limit describes the Darboux III space with parameters 
c = gfr, i) = ^jT, A = 0, e — ;■ — e. For each family of M.S. spaces the second step consists in seeing 
if the Hamiltonians [2] and [3] defined on the Darboux spaces admit a potential V which keeps 
the quadratic superintegrability of the system. An answer to this question was provided more 
recently by [20] [21] who classify the M.S. potentials which we can consider without breaking 
the quadratic superintegrability of the free motion on Darboux spaces. What turned out from 
that classification was that all of the M.S. Darboux systems with potential can be obtained 
through the Coupling Constant Metamorphosis (CCM) applied to M.S. systems on spaces of 



constant curvature. To be selfcontained let us recall briefly what the CCM is. For our purpose 
the coupling constant metamorphosis can be briefly summarized as follows. Let us consider T as 
the kinetical energy term given by T = Y2i j 9^'^PiPj ^-iid U , V are potentials terms independent 
of the arbitrary parameter fi. 

H = T + V-fiU. 

Let S = S{fj,) be an integral of the motion. We deflne the action of the coupling constant 
metamorphosis on H and S as : 

H=^iT + V-E),S^SiH), (15) 

then S is an integral of the motion for H. Let us apply this machinery to the generalization of 
the Kepler problem on a space of constant curvature. The Kepler problem on a two dimensional 
space of constant curvature in cartesian and polar coordinates is given by: 

H = {p^+p)-ti +4^5 (16) 



2 ^^" ^ r 

let us consider the Levi Civita transformation 



H = ^^ ^: ' ipl + ^) - ^^'—^ + A^,5 (17) 



X 2 ' 2 



I y = xy; 6 = 26 
The Hamiltonian in the new cartesian coordinates is 



(18) 



((x2 + y2)-l_A2(x2+y2))2(~2^^2) 

^ iPx -^ PyJ 



i^= "^ ' " "^; ""'"-^ ' ^ ^ {pl+pf,)-l.i{x' + f)-' + X\x' + ff-25) (19) 

or in polar coordinates 



H = ^'"' f"'^'"\ p'r + 4) - f^i^'" + ^'-' - 2^) (20) 

2 r^ 

where A^ = -| , // = | . 
Let us apply the coupling constant metamorphosis to this system. The potential U turns out 
to be 

U = (r-2 + aV^ - 25) 

and the Hamiltonian H is 

fr (r-2-AV)V^ o pj E 

2(r-2 + A2r2-25)^^'^^r2^^r-2 + A2r2-2,5' ^ ' 

It is straightforward to see that the kinetical part of the CCM Hamiltonian H turns out to 
describe a particle moving on a Darboux type space [HI when e = 1, moreover the CCM provides 
directly also the family of possible M.S. radial potentials admitting bound states on Darboux 
spaces in full agreement with the classiflcation given in [20j [21j . 



3. Beyond the quadratic M.S. and Bertrand systems 

In the previous section we have introduced a class of 2-diniensional M.S. systems admitting three 
constants of the motion at most quadratic in the momentum in Euclidean and Non Euclidean 
spaces. As previously said the M.S. entails, for a classical Hamiltonian system, that all of its 
bounded trajectories are periodic, and this is indeed the case for the systems ()171 \2T\i which 
describe elliptical trajectories. Reversing the problem we can look for superintegrability among 
systems admitting stable periodic trajectories. The idea of classifying Hamiltonian systems for 
which all bounded trajectories are closed dates back to Bertrand theorem [IJ proving that the 
only radial potentials with this property are the harmonic oscillator and the Kepler potentials. 
That theorem has been recently generalized to non Euclidean radial spaces by Perlick in his 
remarkable paper ^. Analogously to the original Bertrand's theorem, Perlick found only two 
couples of (metric spaces, potentials ) admitting periodic trajectories for any bounded motion, 
which can be considered as the non-Euclidean deformation of the Kepler and Harmonic oscillator 
on a flat space. Let us report these two families in a conformally flat reference frame as earlier 
obtained in 1131. 



d^' = ;j^-W^k^(dr'+r'd9'+r' sin' edcp'); Vj = -f,{r-^ - kr% /3 G Q (22) 



7.-27 _|_ \2 27 _ 2S n 

d^h = ,2^,-2, _x2r2.)2 idr' + rV + r'sinHd4>'); Vn = ^,-2, + ^^,2, _ 25y ^ ^ ^- (^3) 
The associated Hamiltonian systems turn out to be: 

^2(^-27 _ A2r27)2 ( ^ L'\ f. 



2(r-27 + A2r27 _ 25) V r"^ ) (r-27 + A2r27 _ 25) 

where L2=p2 + _£|_^ ^p^^) = r^.^n^n eN. 
As expected for Hamiltonian systems admitting periodic trajectories the systems (|24l [25]) are 
M.S. with integrals of motion polynomial in the momentum of arbitrary order r = m+n as proven 
in [8], and define implicitly the most general classification of classical M.S. systems with radial 
symmetry admitting periodic trajectories (namely bound motion). This classification of M.S. 
systems overlaps with the previous classfication of quadratic M.S. systems on two dimensional 
spaces whenever we fix the plane of the motion. Let us consider the Hamiltonian systems <\2A\ 
[25]) in the case p^f, = Q. This choice makes the Hamiltonians <\2A\ [25]) comparable with the two 
dimensional Hamiltonians obtained in (jl7|, [2T]) , and in fact they coincide with the (j24l [25]) if we 
consider the particular case /3 = 1, 7 = 1. 
Remark 

The choice of a particular plane of the motion is not restrictive in order to analyze the dynamics 
of a radial system, nevertheless the reduction of a 3-dimensional system to a two-dimensional 
one changes deeply the geometric properties of the space, and in fact the metric space ([22]) has 
a scalar curvature which turns out to be in general non constant : 

Ri = 2{l- p''){r-^ + krP)'' + 2Afk (26) 

while if we consider a fixed S the reduced first fundamental form 



ds^ = o, / , «^o (dr"^ + r^dO^) (27) 

has associated a scalar curvature which is constant for each value of (3 

Rir = 8/3'^k. (28) 

The main consequence of this fact is that the two dimensional reduction of (j24p can be recast 
in ()17p through a canonical change of variables 



r o,-p^ = a—rPr 
^;pe = ape' 



a = /3 = 7 (29) 



H = p^L^!--L^Ll_ (pI, + ^,]- Mr'-i - kr'). (30) 



,r'2(r'-^^_fcr'^2 / ^2 

The same transformation (|29p can be used to recast the (j25p into the (|2ip 

-" ^ 0(^1-2 I \2^/2 _ 9A^ I ^^' "^ r-'2 I + /'^/-2 , \2^/2 _ 9A^ ' ^ '' 



2(r'-2 + A2r'2 - 25) V ^'V (^'"^ + ^^^'^ - 2^^) ' 

This entails that the quadratic M.S. systems on space of constant curvature and its coupling 
constant metamorphosis partners describe the dynamic of the whole family of radial systems 
whose bound trajectories are periodic. 

4. Higher order constants of the motion for Perlick's Hamiltonians 

The above considerations suggest a connection between the quadratic constants of the motion 
of (J17p and (j2ip and the constants of the motion of the Perlick's family which, generally are a 
polynomial of degree n in the momentum variables. Let us analyze this point computing all the 
constants of the motion for the system (J17p . 

4..I. A " Laplace- Rung e-Lenz vector" for the quadratic M.S. systems 

As previously stressed, the system ([T7]) plays a crucial role in the classification of M.S. systems 
with a radial symmetry: in the previous section we have seen how the dynamic of the system 
(|24p reduces to the dynamic of (|17p through a canonical change of variables, moreover the whole 
second family (j25p can be obtained through a coupling constant metamorphosis of (J17p . The aim 
of this section is to analyze the constants of the motion associated to (J17p and to undersatand 
how the transformation (j29p can generate the higher order constants of the motion associated 
to (p^ , [25]) . The two dimensional Kepler system on a space of constant scalar curvature 



J2(^l-l _L u^i\2 / ^2 



H = ''^^'■'~" + ^'^')' ^pI, + ^J _ ^(/-i _ kr') + 4^5, (32) 

is M.S. since it has three constants of the motion: the Hamiltonian itself, the angular 
momentum coming from the radial symmetry pg, and a " Laplace- Runge-Lenz" vector. Let 
us compute the " Laplace-Runge-Lenz" vector from the solution of the orbit equation. 
The orbit equation associated to the system ([32]) is 

de = ^"^"^ L = pg (33) 

this can be recast in the form 



de = K r^ I _ ^3^^ 



(35) 



the equation syinplifies with the change of variables: 



r 



— Ldu , , 

d^ = , (36 

V^2^ + 2/xn - 8/x(5 - L'^u^ - AkL^ 



which can be readily integrated to yield: 



r 2 \—kr^ 

y/2EL^ - 8/i(5L2 - 4kL^ + ^2 

..... (1 + A;r2)Lp, 

sin(6' — t'o) ^ 



V2^L2 _ 8^5^2 _ 4^^4 + ^2 

The orbit periodicity induces the following complex " Laplace-Runge-Lenz" constant of the 
motion 

(1 _ hr^ \ 

l2 fi + i{l + kr^)prL\e-"'^ (38) 

where const = ^JlEL? — 8fi6L'^ — AkL^ + fi"^ is a constant of the motion since it is a function 
oi E,L . Both the Real and Imaginary part of ()38p are constants of the motion quadratic in the 



momentum in radial variable and in cartesian ones, with L = xpy — ypx, Pr = ^'^i' ^^^ '■ 

S = const e^e(r,p.,/f,p,)g-ie ^ f ^2}^^]^ - /i + i(l + kr^)prL] ^^^. (39) 

\ r J r 

4.2. embedding a two-dimensional radial system in a higher dimensional space 
The two dimensional M.S. system ()17p analyzed above can be embedded in a higher dimensional 
space without changing its integrability properties. This operation can be performed if we look 
at the system ()17p and its constants of the motion ([391 , pe) as elements of a universal enveloping 
algebra of a Poisson-Lie algebra whose two dimensional realization can be regarded as just a 
particular choice. Let us introduce the following sl{2) Lie coalgebra 

{J3,J+} = 2J+, {J3,J-} = -2J_, {J_,J+} = 4J3 (40) 

equipped with the trivial coproduct A 

A(1) = 1®1 A{Ji) = J, (g) I + 1(g) Ji i = +,-,3 (41) 

the action of the coproduct defines a homomorphism for the algebra 

{A(J3), A(J+)} = 2A(J+)- {A(J3), A(J„)} = -^(J_); {A(J_), A(J+)} = 4A(J3) (42) 



Let us define moreover A()" as the iteration of the coproduct A()" = A(A(. ...)), then if we 



fix the maximum value of n to n < A'^ we recover a Lie algebra which can be decomposed to a 
direct sum of algebra (P0|) sl{2) (B sl{2).... (B sl{2) whose structure turns out to be the following 

^ V ' 

xTV 

■{A(j3r,A(j+y} = 2A(j+)\ i<j 

{A(J3r,A(J_y} = -2A(J_)\ i<j (43) 

[{A(J_)',A(J+y} = 4A(J3r, i<j 

where i,j < N. From this it is straightforward to verify that the three generators 
A(J+) ,A(J_) ,A(J3) commute with N quadratic Casimirs induced by the action of the 
coproduct on C = J+J_ — J^ 



A(C)' = A(J+)'A(J_)' - (A(J3)')2, 1<1<N. 



(44) 



The elements of the above algebra can be used as fundamental bricks to express the radial 
Hamiltonians and its constants of the motion. Let us introduce a symplectic realization for (|40p : 



D{J_) = x^, D{J^ 



pI^ 



DiJ-s) = xp^, {/, g} = Y^ difdp^g - digdpj 



(45) 



TV 



AN) 



N 



D{j[^^) = Y^xl D{J^:^) = Y^pI, Di4^^) = Y^x,p^^, 



(N). 



N 



(46) 



4 = 1 



4 = 1 



4=1 



We are finally ready to express the system pT|) and its constants of motions as function of 
the generators (jl5|) . Let us put 



f 2 

pi-- 


= C(2) 


r^-. 


-/'^ 


r)_ - 


4^' 



-ie 



Pee = [xpy - ypx)—r- = 
substituting (|^7|) into (fT7|) and ([55|) we obtain 



xVc^ -iixji-"' -p^j):') 



7(2)^ 



H 



(1 + A:ji^))2 (j(2))2 c(2) 



l-kJ 



(2) 



J 



(2) 



J 



(2)- 



M- 



J 



(2) 



+ 4/i(5 



(47) 



(48) 



S 



c(^) ^ ^^ -ti + i{i + k/'^)J^Vc^ 



A 



J 



(2) 



j; 



(2) 



\/c(2y 



The quantity \/A is defined as: 



(49) 



A = Pee' 



{xpy - ypx 



x-iy xVc^ - i{xjf^ - p^Jr) 



r(2). 



J 



(2) 



(50) 



We have used the square root to emphasize that any element of the form ax + f3px can be 
regarded as the square root of a function F{J_ ,J^,J^ ). 



For exampe we have {ax + Ppx)"^ = a'^J_ + 2al3J^ + /3^ J_j_ . The complex constant of the 
motion ()49p provides two constants of the motion given by the real and the imaginary part : 



M.) - ^ {^^^l±^p^ - (£l!Hi^ - ^) ..f - ..i=.,) ,., 



(2)., (2)^ /^(2)n _ 1.7(2) > 



C(2)(l + A;ji'^)J^^^x /C(2)(1-A;J 



^ l(xjf -p.ji')); (53) 



In (j53p we have multiplied the imaginary part by the constant of the motion vC to stress 
that the integral (j49p provides two integrals of the motion polynomial in the momentum. In 
particular the real part turns out to be quadratic in the momentum and in fact is the candidate 
to be the proper Laplace-Runge-Lenz vector. It is possible to verify, by a direct calculation, 
that {iJ(j(2)), (iie(S(j(2), j(i))))2} = by means of the algebra (03]). Since the coproduct A 
defines a homomorphism for the algebra it is immediate to see that 

A^{/7(j(2)),(i?e(5(j(2),j(i))))2} = {F(jW),(/?e(5(jW,jW)))2} = holds for any 
dimension N. More explicitly we obtain that: 



Poisson commutes with: 






t(^)\^. n _ I. TW\f^. tW _ ^ tW\ /W ^^^ 



Li = Re{b)=C' -^^ + -^^ 



j(N) j(N) A-(^ 



(55) 



J'_ 



(1 - k:K.'^)p^xi + 2A;(x • p)^^! - (1 + k:x.'^){x ■ p)px, - ^^^ 



This turns out to be the first component of the Laplace-Runge-Lenz vector of the generalized 
Kepler system on a N-dimensional space of constant curvature [22 , which in the flat case reduces 
to the well known expression for the Laplace-Runge-Lenz vector on E''^ 

£i = p xi - (x • p)px^ - 



Finally it is possible to obtain all the other components of the Laplace-Runge Lenz vector 
through rotation, considering the fact that any radial Hamiltonian Poisson commutes with the 
set of angular momentums Lij = XiPj — XjPi 

Ci = {ai,Li^i} (56) 

which is consistent with the fact that in a A^dimensional space we have A^ possible realization 

for J(i). 



4-3. A "Laplace- Rung e-Lenz" vector for higher order M.S. systems 

In the previous section we have introduced a procedure for the embedding of a two dimensional 
radial system in a higher dimensional space of arbitrary dimension N. The crucial observation 
in order to obtain the constants of the motion for the entire Perlick I family ()24p is that its 
projection on a two dimensional space coincides, up to the change of variables ([29]) . to the 
quadratic M.S. system (fT7|) . This observation can be fruitfully used in the opposite direction, 
namely to obtain first the integrals of the motion of the 2-dimensional version of (|24|) from the 
integrals of p7|) and then to embed them in a higher dimensional space through the application 
of the coproduct A. Let us show how to do this explicitly step by step. First let us apply the 
radial change of variable (|29|) to the system (jlZp and its two dimensionsal " Laplace- Runge-Lenz" . 



(39' 



[P9 = ^Pe' 



(57) 






S = {pI^-^ - /. + i(l + kr')prpe ) e-^' 

P% 1 - kr^ , VI I I, 2n Pe' 

-p2^ ^^ + ^{l + kr )pr — 

which in terms of the coalgebra generators turn out to be 

j{2) ^^j{2)' ^ r^ 

S=\^ ^ ;— -^l + i{l + k/'^) "^^"^ \{^^] . (59) 





Let us remark that, in spite of the fact that the transformation (j57|) is just a change of variable 
for a two dimensional system, this induces a new coalgebraic system (|58|) which is indeed different 
from the original one (|54p . Moreover let us stress that after the angular transformation (|57p the 
constant of motion S is no longer polynomial in the generators J^jJ^, namely the ones containing 
the momentum variables; however /3 = — is a rational number by hypothesis, therefore it is still 
possible to recover a polynomial constant of the motion considering a proper power of S 

S = {Vc(^rS'' = B''(^VAy , m,neN (60) 

B = -^ =^ - A^ + i 1 + k/y)^^^ 



At this stage the constant of the motion S is not yet polynomial in J-(_, J3 since both B and 
A depend on yC, however this dependence is of the type: 

B = {a + ihVC), VA=id + ifVC) 



This entails that if we consider separately the Real part and the Imaginary part of S we get 
two polynomial constants of the motion in the momentum: 



. + ih^/cY{d + ifVcy 



JRe{S)=:F{J+,J_,J- 



\lm{S) = Vcg{J+,J^,J3) 



{H,T} = {H,g} = 0. (61) 



This provides us the way to obtain the "higher order Laplace- Runge-Lenz vector" in any 
N-dimensional space. In particular let us consider the three dimensional version of 



H 



[l + kJ^iJ^^ C(3) 



1-kX 



(3) 



J 



(3) 



/32 jI^) ■ 



^- 



J 



(3) 



AfiS 



which in spherical coordinates : 



xi = r cos a cos ( 
X2 = r cos 9 sin c; 
X3 = r sin 9 



(3) 

(3) 

+ 

(3) 



C(3) = jf j(3)_(j(3))2^p2^ 



xf + X2 + x^ = r^ 

2 2 

Pxi + Px2 +Px3 =Pr + ^ + ^2 sin^ , 
XlPx-, + X2Px2 + X3Px3 = rpr 
o2 



sm 



takes the form 



(62) 



(63) 



H 



[1 + kr^) 



2\2 



Pr + 



1 



ipl + 



:) 



1 — /cr^ 



fJ'- 



2 y ' j3^r^ - - sm^H' I r 

This coincides with (j24p after applying the radial change of variable 



+ 4fx6. 



(64) 



1 ( r'^r'-^ + kr'^f 



1 



Pr' + ^{pI + :tJ-,) 1 - ^'(/-^ - kr'P) + 4/.'5 1 (65) 



sm 



where // = ■^. 

This fully agrees with the results obtained in [^ in which was explicitly proven the M.S. of 
the three dimensional systems (|24|) and (|25|). 



5. Higher order M.S. systems beyond the radial symmetry 

Up to this point we have shown how the entire Perlick 1 family can be obtained from the 
more abstract system ()62p by considering the representation (|46p . However the ()46p is not 
the only possible one, in particular it is possible to get new M.S. systems by considering new 
representations whose symmetry is different from the radial one considered so far. Let us choose 
the two dimensional representation of the algebra sl2 ()46p 

(66) 



D{J'_^>) = x\^xi, D{JX') = vl, +PI2, D{jP) = xwx, + X2Vx,. 



Let us now recast such a representation in radial variables: 

\xi = r cos e, p^^ = cos Bpr - ^pe 
\x2 = r sin 9, p^^ = sin Opr + ^^Pe 



(67) 



r 



^2 



D{J^:^^) = r\ D{j';^)=p'^ + ^, D{jP)=rpr, D{C^)) = pi (68) 

The Casimir representation coincides with the square of the angular momentum pe, and 
the variable 9 can be considered as an ignorable variable since for each algebra element holds 
dgD{J^'^') = 0. As a consequence we can reduce the dimension of D{J^'^') considering pg — )• 6 as 
a fixed parameter 

5(J«)=x2, D(jf)=p2 + ^^ ^(j(l))^^p^ (69) 

which can be considered as a generalization of the representation ()46p when 6 7^ . Of course 
it is still possible to generate the two dimensional representation applying the coproduct A 

A(D(J«))=x? + xi, A(D(J«))=p2^+p2^ + 4 + 4' A(^(J«))=xip.,+X2P.,. (70) 

Let us plug in such a representation in the system ([58|) 

^ (l + fcji^))%(jf)2 C(2) ^ l-kJ^^ ^ ^ 

{l + kr^Y f 2 Pe bl bl \ 1 - kr^ ^ ^ 



The presence of the new terms whose coupling constant is bi breaks the radial symmetry of 
the previous Hamiltonian obtained through the representation ()46p . If we reintroduce again the 
angle 9' = (39 we obtain explicitly the TTW system on a space of constant curvature, which 
from this perspective can be regarded as the two dimensional reduction of a radial 4-dimensional 
system: 



(1 + kr^f ( 2 Pi_\ {l + kr^f ( bj bl \ 1 - kr^ 

2 \Pr+ ^2)+ 2 V/?V2cos2/3-ie'^/3V2sin2/3-i0', 



(72) 



which indeed is M.S. as proven in [23) 

6. The quantum version of Perlick systems 

The analysis performed on the classical radial superintegrable systeins (j24p . (j25p shows how 
the Hamiltonian and its constants of the motion can be viewed as functions of elements of 
the algebra (03]). Its symplectic realization (06) , [70P can provide a wider family of classical 
superintegrable systems. The goal of this section is to take advantage of the above analysis to 
define a quantum realization of this algebra with the aim of inducing, through the coalgebra 
approach, a superintegrable quantization for the systems (pip . (f25|) . As in the classical case, 
the first step in order to realize the above project, is to obtain a two dimensional quantum 
superintegrable version of ()17p . from which, we can obtain the quantum version of ()58p . 



To begin with let us consider again the 2dini Kepler system ()17p . In order to obtain its quantum 
analog we have to deal with the unavoidable ordering problem coming from the canonical 
quantization process that always arises in a non-flat space. For the Kepler system on a space 
of constant curvature the order ambiguity can be solved by introducing the Laplace Beltrami 
operator (the covariant version of the standard Laplace operator) as the quantum counterpart 
of the classical kinetic energy term 

g'^'p^PJ ^ =^di{^9'''d,). (73) 

V9 

In the case of a conformally flat two-dimensional system ds'^ = f{r){dr'^ + r'^dO'^) it takes a 
particularly easy form which, for the sake of coherence with earlier papers, we define to be 
the "direct quantization", namely the quantization obtained prescribing a Hamiltonian operator 
selfadjoint in the Hilbert space L^(M''^, /(|x|)dx) 

1 ^,J,4U^V-^(.J + 1*.^V ,74, 



/(r) \^ r"^ J f{r) f{r) \ ^ r r^ 

The above considerations leads to the definition of the following quantum version of ([T 



H^.n^il±fl(a;^U^m-,'-^. ,75, 

2 \ r r^ J r 

The system ([75]) turns out to be exactly solvable (we can calculate the spectrum and 
eigenfunctions) and above all it is M.S., namely it commutes with two differential operators: 
the angular momentum —ihde and a second differential operator which can be considered as 
the two-dimensional quantum version of the Laplace-Runge-Lenz vector. In the clssical case, 
we obtained the Laplace-Runge-Lenz vector from the solution of the orbit equation ()37p . In 
the quantum case the exact solution of the bounded spectrum and its eigenfunctions can be 
computed algebraically factorizing the Hamiltonian operator in terms of two "ladder operators" 
H = a) a . On the same footing of the classical case the quantum Laplace-Runge-Lenz can be 
obtained by means of the same tools already used to working out its exact solution, namely 
using the properties of the ladder operators. The idea of using the ladder operators to get the 
extra integrals of the motion for a superintegrable systems is not new see i.e. [21] j25j . In the 
following we will adopt the same philosophy to compute the quadratical constants of the motion 
for the two-dimensional system with the perspective of generalizing this to the N-dimensional 
case as already seen in section Q for the classical case. As a first step let us separate variables 
in the Schroedinger equation, i.e. putting ij) = e^^^ p{r) 

Hip{r) = [-n^^ll^ {dl + \dr - ^) - /.i^) P{r) = mr). (76) 

The radial Schroedinger operator ([75]) can be factorized introducing the operators a|,d/ : 

(77) 









Hi-Ei = a]ai (78) 



where Ei = -^^^ + 2kh\l +'^f-^. 



Let us stress that the operators introduced in (f77|) satisfy the following important relation: 

aia\ + Ei = Hi+i = a\^-^^ai+i + Ei+i. (79) 

Such a property is known as shape invariance and it is an integrability condition ( for a review 
on SUSY QM see i.e. [2^). It provides a way to obtain the entire discrete spectrum and the 
eigenfunctions. Let us show briefly how the shape invariance condition works. 
Let us consider the operator aj, ^a^+i, and the radial wavefunction po,/+i which is annihilated 
by ai+i ■ 

the wavefunction /Oo,i+i can be regarded as the ground state of the Hamiltonian 

[ili+\ - -Ez+Opo.z+i = al+iai+i/Jcz+i = 0. (81) 

Using (j79p equation (|8ip can be recast in the form 

Oiolpci+i = [Eij^x - Ei)po^i+i. (82) 

Multiplying from the left by a| provides the first excited state of the Hamiltonian a/iai 

a\ai{a\pQ,i+i) = {Ei+i - Ei){alpQ^i+i) (83) 

or equivalently 

Hi{alpo,i+i) = Ei+i{ajpo^i+i). (84) 

It is then straightforward to obtain all the bound states and their energy eigenvalues by 
iteration 

n-l 
Pn,l OC JJ al_^_jPo^l+n] En,l = En+l- (85) 

i=o 
This provides the solution of the two dimensional eigenvalue problem ()75p 






H^l^^^i = He'''pn,i = -^^TTTT^ n? + 'Ik^'il + ^ + o) - ^^ V-n,;. (86) 



Let us remark that the spectrum shows the so called "accidental degeneracy" namely it 
can be described by the single quantum number M = n + I. This fact can be regarded as a 
direct consequence of the shape invariance property ([79]) . The operators ai and aj act as ladder 
operators and connect different " iso-energetic" radial eigenfunctions : 

{aiPn,l fX Pn-l,l+l .„„x 

alpn,l+l OC Pn+l,l- 

The radial ladder operators ([57]) . as introduced above, are defined for a specific value of I. 
We can generalize their action on a generic eigenfunction ipn,i by replacing the quantum number 
/ by the angular momentum —ihdg. 



Let us define the operator S: 

S = ^ (-ih{l + kr^)dr + '-^^ - h^—^de] . 

It is now just a matter of computation to see that S acts as a ladder operator for the set of 
eigenfunctions il)n^i '■ 



—^ -ih(\ + kr^)dr H 

V2 \ —ihde 2 



hnx^i = ^ -^Kl + kr')dr + '^ , - h^-J^de e*('+^)^p. 



it is straightforward to see that its Hermitian conjugate S^ behaves as a : 

5t = ^ (-ih{l + kr^)dr - . '^ , + h^^^de] (90) 



\/2 \ -i^iae + ^ 



2 



r 



^/^ , fel-^'^^^ \ jw. 

71 r^'^'^"'^' ^^'' " ^^^^5^T^ 



5t^„,, = ^ -^^(1 + fer2)9, - '!" . + ^i^^^^Se e^'Vn,/ 



i(/+l)6» y^ ^ 1 _ jt^2\ 

ih{l + /i;r^)9r - — y\ + ^^^ /'".^ °^ V'n.-i,«+i- 



^/2 \^ ' ' /i(Z + i) 

On the basis of the above considerations, it is natural to define the operator 

5 = iS{-ihde - ^) (91) 



which can be regarded as the quantum version of (|49p . Furthermore the definition given in 
(|9ip provides a straightforward way to obtain the two dimensional quantum Laplace-Runge-Lenz 
vector 

^ = -^^ (92) 

whose classical limit coincides with the ([?T|) 

6.1. embedding a two-dimensional radial system in a higher dimensional space 
Given the quantum version of (J17p and its Laplace-Runge-Lenz vector ()5ip the second step 
consists in describing them in terms of the algebra (j43p . In order to accomplish to this task let 
us introduce the following quantum representation of ([I3|) 



[Js, J+] = 2ihJ+, [Js, J_] = -2i/iJ_, [J-Js] = 4ihJ3 (93) 

[A( J3)\ A{J+y] = 2ihA{j\)\ i < j 

[A(J3)S A(J_)^] = -2inA(J_)', i < j (94) 

[[A{j^y,A{j+y] = mA{jsy, i<j. 



This time the Casimir operator is the symmetrized version of the 



A(C) 



J, _ A(J+)^A(j5) + A(J_)^A(J+)^ 



(A(J3 



\N\2 



A quantum reahzation of ()94p is given by the foUowing differential operators: 



N 



i{N) 



N 



(N), 



N 



D{J^_^^) = Y^xl D{ji^^) = -n'^dl, D{J^/^) = ^{-^hx,^., 



Nih 



1=1 1=1 

Considering the quantum version of ()47p 



1=1 



r^ = jl 



-ihdr = Pr 



1 f I ih 



/J(2) 



(2) 



e pg 



f{2) 



V^ - ^^(xjf ^ - f^p^ + ihx) 



it is immediate to obtain the "algebraic" version of ([75 



(1 + yfcJi^ 



H = -fv 

(2)^2 



2(1 + A:r 



2\2 



1„ <9l\ l-kr'^ 



d., + -dr + Hj I - ^- 



J 



1 ' 7(2)^2 , 2i?l j(2) ^ _J_q{2) 



(2) 



{jrr + 



1-A;J: 



(2) 



/ 



(2) "3 



J 



(2)' 



A^- 



J 



(2) 



(1 + fc/ 



(2)^2 



i?^-^- 



i-fcj' 



(2) 



J 



(2) 



and of its constant of the motion 



(95) 



(96) 



(97) 



(98) 



C 



S + S^ 



i-kfy)j^:hmkji^^+2HjPf-^^ 



Ci = xi \-kh^ + 



(99) 

The operator £» can be regarded as the square root of a function J^(J+, J-, J3) in the sense 
that : 

Cj = J-(J+, J_, J3) : [-F(J+, J_, J3), iJ] = 



which in view of (|94|) commutes with if, as previously stressed for its classical version 



It is also possible to evaluate directly the expression [£«, H] if we enlarge the algebra (jMj) with 
the two elements of the Heisenberg algebra px^ , Xi : 



[Px,,J. 



j{N) 



0, 



Xi,J\_ '] = 2ihpx 



[Px^,J3 ] = -ihPxi, [xi,J\. '] = -ihx 



(100) 



[Px,,J. 



{Nh 



^ Li L<Xj • iXj 7 1 " __ 



Wi 



0. 



This makes it straightforward to obtain the N-dimensional version of (|75p by choosing the A^ 
dimensional realization of the algebra generators : 



H 



(^) _ (l + fcjj^^)^ -(^) _ ,. 1-fcj! ^) 



-JX'-^^- 



J 



(N) 



MN)__ 



(101) 



(102) 



Xi \-kh'^ + {l-kf^^)jf''+3ihkJ^f'^ + 2k{J^ 



fW^2 ^ 



J 



(N) 



[h(''\c 



Wi 



0. 



l + A:j!^Vr + f-^i^^')p, 



(103) 



6.2. Higher order quantum "Laplace- Rung e-Lenz" vector 

Carrying on the parallelism with the classical case we can get the quantization of ()58p from that 

of (HZD applying the substitution 6 = PO' to ([75]). 



H = -h- 



H' = -K' 



2(1 + A;r 



2\2 



dl + -dr + 



dl 



l-kr^ 






2n2 



(l + kr^) 



dt + -dr + 



/32r2 



1 — fcr^ 



A*- 



(104) 



S 



le 



Ad 



Pi 1 j^T h 

, ih{-ihde - -)(1 + kr'^)dr + ifj,- h {-ihdg - -)dt 

V2 V 2 r 2 



5' 



ze 



^ ^ ^ni^^e, - ^)(1 + A:r2)5. + ^/. - ^^-^i^de^ - ^)5,, ) . (105) 



As in the classical case, this angular substitution doesn't break the M.S. whenever we consider 
/3 G Q. This fact becomes evident if we take in account how the angular substitution acts on 
the spectral properties of the new Hamiltonian . Let us compute the new eigenfunction: 






r 

1 



E{n + l) = E{n + ^) = E{n + ^^^^), /3 = — , mi, ma G N. 
p m2 m-i 



(106) 
(107) 



The new bound spectrum ()107p still exhibits the "accidental degeneracy" having an infinite set 
of isoenergetic eigenfunctions 

< V'njI-^VnJ > = < ll^'n+smul-sm2\H'Wn+3m^,l-sm2 >' ^^ ^ I^- 

The operator S as defined in ()105p is no longer a "ladder operator" since it loses the property 
of connecting isoenergetic eigenfunctions of (|75p . however it is possible to restore such a property 
redefining a new operator ^S* as an appropriate power of S : 



Ap-im2e "^1 / 1 _ t, 2 _ -t \ 

jS = 5™i = — ^ n ( -^^^i(l + kr^)dr + ¥- i^-^V,{—de - hj)) , (108) 



The expression (jlOSp can finally be used to define the additional constant of the motion ^C 
for (flM)) 

.£ = ^^^. (109) 

The crucial point is that the Hamiltonian operator ()104p if expressed by means of the relations 
([97|l turns out to define a new coalgebra operator: 

2 \^j(2)' 3 y ^^2) ^2j(_2) ^2/2) y ^ /J(^ 



As we will prove in the following sections, the operator (jllOp embedded in a A^ dimensional 
space, turns out to be exactly solvable and its bounded spectrum exhibits the so called maximal 
degeneracy characteristic of the M.S. quantum systems. The proof of the M.S. for the system 
(jllOp in any dimension N requires, from a coalgebraic perspective, the evaluation of the 

commutator [H'^'^\ pt {{jT''}:{j'f'}^Px,x\] = in terms of the algebra ( [Ml , [10(1 . Such 

a computation is in general very hard and needs, for higher order constants of the motion 
(/? 7^ 1), several days of machine-time, however we have verified this explicitly in the case 
(/c / 0, /3 = 1; /c 7^ 0, /? = |; A; = 0, /3 = ^, A; = 0, /3 = 2) namely up to the third order constants 
of the motion. If we conjecture that the M.S. holds for any value of the rational parameter (3 
than it is straightforward to obtain the explicit expression of the constants of the motion in any 
dimension N . Let us prove this considering for simplicity the square of the Laplace-Runge-Lenz 

f^C = Q i{Ji },{Jj }), which in analogy with the classical case can be expressed by using 



only the generators of ([9 

[^(U^'^}),^({4'^},{-/r})]=0^[i?({j(^)}),a({j(^)},{j("^)})]=0, m<N. (Ill) 



Let us remark that if we set m = 1 then there are A^ possible reahzation for the set of elements 
{Jj } namely 



jw = _h2Q2^^ j|ij ^ _inxid,^ - -, A'> =xl 1 = 1, 2, ....N 



providing N constants of the motion {Q{{J^^'^t{J^^'^i)) which, together with the set of 
angular momenta, makes the quantum system (jl04p a M.S. Hamiltonian. 



6.3. spectrum and eigenf unctions of the n- dimensional system ^104 ) 



Let us conclude the analysis on the N-dimensional embedding of ()104p evaluating explicitly its 
eigenfunctions and the spectrum 






A;J^ 



(112) 



let us introduce the N-dimensional radial coordinates : 



Xj = rcos6j TT sin^fc; xn = r Y\ sin^^ (113) 

k=l k=l 



consequently the generators ([M]) become 



^2 ' 3 ''2 



jW = .^ jW = -;,^g^2_ »i^v^-^; g^^i_^ jW = _^;,,a.-^ (114) 



Af-1 /i-1 



Using the radial representation ()114p the Hamiltonian operator ()116p takes the form 

(115) 
Let us separate variables in the eigenfunction ip{r,9)n,i,ii...ij^_2 = p{^)n,iY{9)i,ii...iN-2- 
The functions 5^(^)/,Zi..jjv-2 ^^^ ^1^^ ^^1 °1 hyper spherical harmonics functions, namely the 
eigenfunctions of L^ which satisfy the eigenvalue equation given by: 

L^Y{0\i,..,i^^, = hH{l + N- 2)Y[e\i,..,,_,. 
The above factorization leads to the radial Hamiltonian operator: 

'■2^2 / AT 1 f{2in^N_2) h^ ( I \^^^ ^A l-kr"^ 



- ,„, . (1±|^ (_,.a. _ ,.^a, . ^^ffi±^ - ^ (^ - ') '" - 



however the operator Hr is gauge equivalent to ([75 



(116) 






considering the expressions (j85p and (j86p we can finally determine the spectrum of (J116p 

^n i = f -, ^^-^, r- + '^khH,^ + n+ ^^ + -f - — ^ . (118) 

"' ^ 2^2(1 + n + ^ + i)2^ V 2/3 ^2^ 2^ ^ > 

Before concluding the section let us also remark that the Hamiltonian (J116p is indeed 
equivalent to the direct quantization of the N-dimensional Perlick I system defined in 
(|74p . If we define the new radial variable r = r'" and apply the gauge transformation 

l2-N){l-fJ) y,j^. (iV-2)(l-g) , . , TT -T ■ 

r 2 j^K^iJj- 2 we obtam the new Hamiltonian operator: 

r' 2 m^'r' 2 = -h^ — ^ -L ^v^ - /i(r'~'' - kr'P). (119) 

Zp 

This means that the coalgebra symmetry induces the " direct quantization" and for this reason 
the direct quantization is the quantization which keeps the M.S. propety for the Perlick systems 
as already stressed in [27j, [TlJ. On the other hand if we demand a covariant quantization we 
have to consider a quantum correction term to the potential proportional to the scalar curvature 
of the system which makes the Laplace Beltrami quantization of the N-dimensional system ()24p 
gauge equivalent to (|119p and therefore M.S. . 

fi2 1 N — 9 

Hl.b. = - Y^^^ {Vd9'''d,) + F(r) + /i2____i?(^) (120) 



2-JV ~ N-2 



where / — ^ 



f^Hdf^ =Hl.b. (121) 

r2(r-''+fcr'5)2 



7. 4-dimensional Perlick quantum systems and TTW systems 

We conclude the paper by obtaining explicitly the generalization of the quantum TTW system on 
a space of constant scalar curvature by reduction of the 4-dimensional Perlick quantum system: 

H"> = "^"'"'' (^0<y + ^ jf + ^d") + ^^" -P) - ,i^. (122) 



Following the same strategy shown in section 5 let us compute the reduction of the 4- 
dimensional realization of (1961) . 



j^^) = x\ + xl + xl + xl; (123) 

j^^> = -ih{xidx^ + X2dx2 + xsdxs + Xidx^ + 2); (124) 

jf = -h\dl + dl + dl^ + dl). (125) 



This representation can be reduced to a lower dimensional one by using a bi-polar coordinates 
system: 

xi = ricos</>i ; X2 = risin(/)i (126) 

X3 = ^2 cos 02 ; X4 = r2 sin (/<2 



J^^=rl+rh (127) 

jf ^ = -ihindr, + r2dr^ + 2); (128) 

J|^) = -h\d'^^ + I5., + ^,dl + 9,2, + I5., + la^J. (129) 

'1 '1 '2 '2 

Since the generators are independent on the new angular variables 0i , (j)2 it is possible to 
get rid of the two degrees of freedom coming from (pi , (f)2 and to obtain a new two dimensional 
system which inherits the properties of the original 4-dimensional one. Let us show in a few 
algebraical steps how to perform the reduction for the quantum case 

<1p{xi,X2,X3,X4)\H^d\'4'{xi,X2,X3,X4) >= (130) 

flp{xi,X2,X3, X4,)*H4d'l(j{xi,X2,X3,X4)dxidx2dx3dx4 = 

frir2'ip{ri,r2, (pi, (P2)* HM'4'{ri,r2, (pi, (p2)dridr2d(pid(l)2 

the next step is to separate the wave function y/rir2'4){ri,r2,(pi,4>2) = ip{ri,r2)e'^^^'^^ e^^'^'l''^ so 
that the ()130p turns into: 



/ 



f'ip{ri,r2)* H2dip{ri,r2)dridr2 (131) 

where the reduced operator is defined by: 

H2d = Vn^ {J^ j e-''^^'e-''^^^H4de'''^^e''^^^d(pid^^ -L=. (132) 

The reduced Hamiltonian can be now written in terms of the reduced version of the original 
generators 

7(4) , f(4) _ /Z—-(_}_ [ ^~ih<j>l,-il2<t>2li^)jh<t>lA<t>2, 1 



Jl^*' -^ J^> = Vn^(T^ / e~''^'^^e-''^'^^Jr'e"-"^^e''^^^d(Pid(P2) 



4vr2 J * ^Jrlr2 

p)=rl + 4^f^^; (133) 

Jt^ = -ihindr, + r2dr, + 1) = JP; (134) 

-d) 9.^9 1-4/? o 1-4/^, ^(2) ^94/? -1 ^o4/3-l ^ ^ 

jf = -h\d'^^ + -^ + d'^^ + -^) = Jf + h'^^ + h'^^. (135) 

After the reduction, the 4-dimensional representation coincides with the 2-dimensional ones 
except for the generator J+ wich has an additional centrifugal term depending on the quantum 



numbers li,l2 coining from the degrees of freedom we had cut off previously. Eq. ()133p can 
be regarded as a quantum non radial representation of the algebra ()94p . If we plug this 
representation into the "algebraic" Perlick Hamiltonian (|116p we obtain the system : 

fl'^' = ^i^ (-n'ai - n-'-B. - s^l* + -^ + -^\ - ,'-^^ (i36) 

2 y r r^ r^ cos^ |- r^ sm"' j J r 

where we have expressed the generators ()133p in terms of the variables 

9' 

ri = r cos — 

. 9' 
r2 = r sm — . 

We have reabsorbed the "old" quantum numbers "/i" , "/2" in the parameters 6i = . „^^ , 

1-4/2 

^2 = 4gi^ ■ This is a generalization of the the TTW system to spaces of constant scalar curvature. 

8. Conclusions 

The main results of this paper can be summed up in the following three points : 

1) We have introduced a technique based on the coalgebra to generate M.S. classical and 
quantum systems in N-dimensions from 2-dimenionl ones. In particular we have shown how 
the family of 2-dimensional radial M.S. systems with quadratic constants of the motion can 
generate superintegrable systems with higher order constants of the motion of arbitrary order 
by immersion in spaces of dimension N > 2. As a particular case we have obtained the 3- 
dimensional Perlick classification of classical Bertrand-Hamiltonians by embedding of the two 
dimensional systems Darboux III and Darboux IV. 

2) All the systems analyzed in the present paper have been considered both in their classical and 
quantum mechanical version. In particular we stress the connection with the SUSY Quantum 
Mechanics. The radial Schroedinger equation associated to this class of superintegrable systems 
can be exactly solved by mean of the shape invariance technique and moreover all the integrals 
of the motion can be straightforwardly obtained through the embedding of the ladder operators 
in a higher dimensional space. 

3) The connection with the family of TTW systems is also discussed. We show that the family 
of radial M.S. systems and the TTW share the same coalgebra and indeed the TTW can be 
regarded as a dimensional reduction of a 4-dimensional M.S. radial system inheriting in this way 
all its integrability and solvability properties. 
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